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I. INTRODUCTION 

Optically trapped neutral atoms are being actively de¬ 
veloped for cmantum simulation and quantum computing 
applications [ll, 0 and there has been substantial recent 
progress in improving the fidelity of one- and two-qubit 
gate operations [sl-Q- Several different approaches are 
possible for encoding qubits in neutral atoms. For exam¬ 
ple collective encoding provides a method for establishing 
a multi-qubit renter in the collective states of a single 
atomic ensemble [3- One of the challenges in implemen¬ 
tation of collective encoding is measuring the state of a 
single qubit without disturbing the rest of the register. 
This can in principle be done by state selective excita¬ 
tion to a Rydberg level followed by ionization. This has 
the drawback of suffering from less than unity quantum 
efficiency of practical ion detectors, plus the problem of 
atom loss. After each measurement of a bit value of |1) 
an atom is lost and has to be replaced from the collective 
reservoir state. The number of measurements which can 
be made before the reservoir is depleted is thus limited 
by the number of atoms in the ensemble. An alternative 
is to perform a Rydberg gate between the register to be 
measured and an auxiliary register (or single qubit) in a 
neighboring trap. The state of the auxiliary bit can then 
be measured without atom loss. This has the drawback 
of requiring a longer range gate to be performed. For 
qubits encoded in a single atom, optical trap arrays can 
be used to define a multi-qubit register 0, H, [l3- Also 
in this case measurement of the state of a single qubit 
without disturbance of proximal qubit locations is chal¬ 
lenging due to the isotropic distribution of light scattered 
during a measurement. 

State measurements may also be based on cross en¬ 
tanglement of two different atomic elements located in 
the same trap, or nearby traps. By creating entangle¬ 
ment between qubits encoded in different types of atoms 
the quantum state of a qubit encoded in atom a can be 
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FIG. 1. (color online) The resonance condition for interspecies 
dipole-dipole coupling between initial states a, b and target 
states a, l3 is S^a = d/Bb- 


measured via light scattering from the qubit encoded in 
atom b. This is analogous to the mixed species quantum 
logic spectroscopy previously demonstrated with trapped 
ions[ll|. For example Rb atoms have D1 and D2 res¬ 
onance lines at 795 and 780 nm while Cs atoms have 
D1 and D2 lines at 894 and 852 nm. The large separa¬ 
tion implies that measurements, as well as optical pump¬ 
ing and state preparation, can be performed indepen¬ 
dently on nearby atoms. Atoms of different species a 
and b can have a strong dipole-dipole interaction due 
to a Forster type mechanism when the energy defect 
hS = h{Saa + = {Ua “ Ua) + (t//3 “ Ub) is Small 

as shown in Fig. [TJ Here a, b denote initial quantum 
states and a, /3 the dipole coupled states. In this paper 
we provide a detailed analysis of interspecies Forster res¬ 
onances for Rb and Cs atoms and analyze the application 
of the interspecies coupling to quantum non demolition 
(QND) state measurements. 

The structure of the paper is as follows. In Sec. [TTl 
we provide general formulae for calculating interspecies 
dipole-dipole interactions. The formalism generalizes the 
results of [a to the situation where the laser excited 
atoms are not in the same quantum state. In Sec. [1111 
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we present a list of useful Forster resonances for Rb-Cs 
coupling. In Sec. IIVI we list the strongest resonances 
for coupling Rb-Rb and Cs-Cs using different principal 
quantum numbers n for each atom. In Sec. [V]the angu¬ 
lar variation of the interaction is calculated for isotropic, 
and strongly anisotropic cases, and in Sec. [Vl] we dis¬ 
cuss the problem of qubit measurement and show how 
the interspecies coupling can be used for fast measure¬ 
ments with very low crosstalk. Section IVIII summarizes 
our results. 


II. DIPOLE-DIPOLE RYDBERG INTERACTION 
BETWEEN DISTINGUISHABLE ATOMS 


In this section we provide explicit expressions for calcu¬ 
lating the interspecies dipole-dipole interaction between 
atoms a and b leading to Rydberg blockade. Our nota¬ 
tion mostly follows the theory of [l^ with some modifi¬ 
cations, and slightly generalized to allow for the initial 
Rydberg pair states to be distinguishable. We character¬ 
ize the strength of the interaction for a particular angular 
momentum channel by the Cs and van der Waals coeffi¬ 
cients. The label 7 ^ = (za^riaja^ja) denotes the quan¬ 
tum numbers of a single Rydberg level a. The coupling 
{ab) specifies an interaction channel k coupling 

a pair of atoms in fine structure levels a^b to a pair of 
atoms in fine structure levels a^(3. Here z specifies the 
atomic species, n is the principal quantum number, I is 
the orbital angular momentum, and j is the total elec¬ 
tronic angular momentum of a fine structure state. We 
assume single electron s = 1/2 atoms throughout. 

We define the coefficient of channel k as 


C3,k{a,b,a,l3) = 


{la\\ra\\la){lf)\\n\\lb) 
\/ {“^ja + l)(2i/3 + 1) 


( 1 ) 


with = e^/ 47 reo, e is the electronic charge, cq is the 
permittivity of free space, and ( 7 a|ka|| 7 a) is a reduced 
matrix element in the fine structure basis. This differs 
from the notation of [T 2 | where the C 3 coefficient was 
defined in terms of radial matrix elements in the n, I ba¬ 
sis. Note that C^^k depends on a total of 14 parameters: 
Za 1 Zhi •> 3a i jb-) '^a ? i ja •> ^/3 •> ^j3 •> j{3 • 

The energy defect for channel k is h6k = ^(^aa + ^/36) = 
[U ( 7 a) - t^( 7 a)] + [U ( 7 / 3 ) - ^( 75 )] • In fhe approximation 
that a single channel dominates the interaction the en¬ 
ergy shift of a Forster eigenstate \uki) depends on the 
interatomic separation R as 


1 /2~ 

4:Dke{ma,mb)C^,^\ 

^ ^ ) 

( 2 ) 

The angular factor Dki is always positive so for Sk > 0{< 
0) the interaction is attractive (repulsive). The long range 
van der Waals interaction for eigenstate i in channel k is 

jj 1 


TT f \ 

Uk£{ma,mb) = — 


1 - 


We define a crossover distance Rc marking the boundary 
between a 1/R^ resonant interaction and a 1 van der 
Waals interaction by 


Rc 


V m ) 


The angular factor Dk{ma,mi)) depends on the quantum 
numbers of the interacting states and is calculated with 
the method described in Appendix A. 

When we consider the interaction of atoms of different 
types, either two different atomic elements, or two differ¬ 
ent isotopes of one element, we have 7 ^ 7 ^ 75 and only 
include the coupling (ab) {a/3). Also for atoms of the 
same type but with 7 ^ 7 ^ 75 there will usually only be a 
single coupling (ab) {al3) which is dominant. The Dki 
values for channel k and eigenvector i for nsi /2 states are 
given in Table H Interaction of atoms of the same type 
which are prepared in the same levels, 7 ^ = 76 , will have 
two sets of couplings of the same strength: {ab) {a^) 
and {ab) {I3a). This gives the twice larger values 
given in Table U which are in agreement with the values 
given in Table I of 

Starting with a specific molecular Rydberg state j?/^) = 
Cij\mai,mbj) the interaction energy due to channel 
k is found by decomposing into the Forster eigenstates 
\uk£). Writing \ 2 p) = Y.^Ck£\uk£) with Ck£ = we 

have 


^ \Ck£\^Uk£- 
£ 

When there are multiple interaction channels {/c}, cor¬ 
responding to additional values of the situation 

is more complicated and in general has to be treated by 
numerical solution of the eigensystem of the matrix in 
Eq. jAH), extended to include multiple channels. When 
R ^ Rc so the interaction energy is small compared to 
the Forster energy defect there is negligible amplitude 
of the target states a, /3 and in a first approximation we 
may assume the energy shifts are additive. In this van 
der Waals limit the interaction energy is 

^|7>,vdW = ^ \Ck£\‘^Uk£,vdW- (3) 

k,£ 

At small R where the interaction is resonant and there 
is substantial state mixing we must account for coupling 
between channels, which is most conveniently done nu¬ 
merically. The interchannel coupling may lead to nonad¬ 
ditive behavior, as has been discussed previously fl^ 


III. RB-CS FORSTER RESONANCES 

Forster resonances for Rb-Cs coupling occur for a 
range of angular momentum channels. The simplest case 


^ To compare the values in Table I of m with those given here it 
is necessary to account for the different definitions of C 3 . 
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TABLE L Interaction channels k, eigenvalues Dm, and eigenvectors \um) for atoms in nsi /2 states. Eigenvalues Dm are for 
atoms in identical initial states 7a ^ 76 with allowed couplings [ah) ^ Eigenvalues D'm are for atoms in initial states 

7a = 76 with allowed couplings (ab) ^ {olI3) and {ah) ^ (^da) . The eigenvectors specified in terms of states |?Tia,m6) are the 
singlet and triplet states \us) = ^(|l/2, -1/2) - | - 1/2,1/2)), \uto) = ^(|l/2, -1/2) + | - 1/2,1/2)), \ut±) = | ± 1/2, ±1/2). 
States s, to, t± are labeled as i = 1, 2,3,4 respectively. 


channel k 

joe 

i/3 

m = rua -h rub 

\uki ) 

Dki 

D'u 

1 

1/2 

1/2 

0 

|u4 

0 

0 




0 

|wto) 

16/9 

32/9 




±1 

\ut±) 

4/9 

8/9 

2 

1/2 

3/2 

0 

\Us) 

2 

4 




0 

|wto) 

2/9 

4/9 




±1 

|Mt±) 

14/9 

28/9 

3 

3/2 

1/2 

0 

|Ws> 

2 

4 




0 

|wto) 

2/9 

4/9 




±1 

|Mt±> 

14/9 

28/9 

4 

3/2 

3/2 

0 

|Ws) 

2 

4 




0 

|wto) 

34/9 

68/9 




±1 

|Mt±> 

22/9 

44/9 


is excitation of nsi /2 levels which are dipole coupled to 
npi/ 2 ,np 3 / 2 . Figures [21 [3] show the energy defects for 
all four fine structure channels. There are a large num¬ 
ber of resonances with the Rb principal quantum number 
either larger than or smaller than that of Cs. Table [TTl 
lists the strongest resonances for 40 < upb < 90. Radial 
matrix elements were calculated using the WKB approx¬ 
imation ofR5| with quantum defect values taken from 
Refs. [H, fl^for ^^Rb and [H, [HI for Cs. 

The strongest resonance in the table (last row) pro¬ 
vides an interaction strength of 2 MHz at R = 20 jam. 
Even stronger resonances are available at higher n. For 
example the resonance at npb = 121, ncs = 124 gives 
MHz scale interaction strengths at R = 45 jam. Note that 
the energy defect at a resonance can be either positive 
or negative so the interaction can be either attractive or 
repulsive. This behavior is distinct from the intraspecies 
coupling for Rb-Rb or Cs-Cs excited to the same ns states 
for which the interaction is always repulsive (see Fig. m. 





FIG. 2. (color online) Energy defects for interspecies cou¬ 
pling with nscs > nsRb, npRh = nspb, npcs = nscs-1. 
The different curves show the {jpRh, jpCs) channels in the se¬ 
quence (3/2, 3/2), (1/2, 3/2), (3/2,1/2), (1/2,1/2) from top to 
bottom. 


IV. FORSTER RESONANCES OF RB OR CS 
ATOMS 

Analogous to the Rb-Cs Forster resonances studied in 
Sec. [ml there are resonances for Rb-Rb or Cs-Cs in¬ 
teractions. Figure [4] shows the energy defect for exci¬ 
tation of two Rb atoms or two Cs atoms to the same 
nsi /2 state. Even at large n the energy defect is sub¬ 
stantial for the dominant channel which limits the inter¬ 
action strength [HI. The energy defect can be reduced 
using an external field to give a so-called Stark tuned 
Forster resonance, as has been demonstrated experimen¬ 
tally with dc[20| oi* ac[2l| fields. Alternatively, the in¬ 


teraction strength can be increased substantially, with¬ 
out an electric field, by exciting each atom to a differ¬ 
ent n for which there is a Forster resonance as shown 
in Figs. [5l [H This type of resonance has been used to 
advantage in recent atom-photon coupling experiments 
with Rb atoms [HI- Tables [Till [13 list the strongest in- 
traspecies resonances. Comparison of the tables for inter¬ 
species and intraspecies resonances show that they have 
similar strength. 






































4 


TABLE IL Forster resonances for Rb-Cs ns 1/2 states. Only resonances with defects less than 0.0005 times the level spacing 
and |C 3 ,fc| >1.0 GHz/xm^ for 40 < nRb < 90 are listed. The van der Waals interaction strength and Rc distance are listed for 
the strongest eigenvector for each channel, which are |nto) for channels 1,4 and \us) for channels 2,3. 


channel 

Rb|7a) 

Rb 7a) 

Cs|7f,> 

Cs|7/3) 

5/2n (MHz) 

Cs^k (GHz /xm^) 

Uk,vdw (GHz /xm®) 

Rc (/im) 

4 

45si/2 

45p3/2 

47si/2 

46p3/2 

28.0 

-1.29 

-223 

4.47 

4 

46si/2 

46p3/2 

48si/2 

47p3/2 

-11.5 

-1.41 

656 

6.21 

3 

48si/2 

48p3/2 

5Isi/2 

50pi/2 

-5.71 

1.69 

994 

7.47 

1 

59si/2 

58pi/2 

57si/2 

57pi/2 

-16.6 

-3.54 

1350 

6.58 

1 

6Isi/2 

61pi/2 

65si/2 

64pi/2 

2.65 

-4.80 

-15500 

13.4 

2 

68si/2 

67pi/2 

67si/2 

67p3/2 

5.25 

6.5 

-16100 

12.1 

2 

69si/2 

68pi/2 

68si/2 

68p3/2 

-7.40 

6.92 

12900 

11.0 

3 

71si/2 

70p3/2 

69si/2 

69pi/2 

9.35 

8.01 

-13700 

10.7 

3 

72si/2 

71p3/2 

70si/2 

70pi/2 

-7.99 

8.51 

18100 

11.5 

2 

72si/2 

72pi/2 

75si/2 

74p3/2 

4.61 

9.65 

-40400 

14.4 

2 

73si/2 

73pi/2 

76si/2 

75p3/2 

-4.31 

10.2 

48400 

15.0 

3 

77si/2 

77p3/2 

8Isi/2 

80pi/2 

-2.19 

12.3 

138000 

20.0 

1 

81si/2 

80pi/2 

00 

lo 

78pi/2 

6.31 

-13.4 

-50800 

14.2 

1 

82si/2 

81pi/2 

79si/2 

79pi/2 

-6.41 

-14.2 

55600 

14.3 

1 

84si/2 

84pi/2 

89si/2 

88pi/2 

-2.43 

-18.2 

243000 

21.5 


TABLE III. Dominant Forster resonances for Rb-Rb nsi /2 states. The van der Waals interaction strength and Rc distance are 
given for the strongest eigenvector for the listed channel. 


channel 

Rb|7a) 

Rb|7a> 

Rb|76> 

Rb|7/3) 

<5/27r (MHz) 

C3,k (GHz nm^) 

Uk,vdW (GHz jJjJJl ) 

Rc (/im) 

4 

38si/2 

38p3/2 

39si/2 

38p3/2 

4.62 

-0.621 

-315 

6.39 

1 

49si/2 

49pi/2 

51si/2 

50pi/2 

30.3 

-1.91 

-214 

4.38 

1 

50si/2 

50pi/2 

52si/2 

51pi/2 

-31.5 

-2.08 

244 

4.45 

2 

65si/2 

65pi/2 

67si/2 

66p3/2 

32.0 

6.30 

-2480 

6.53 

2 

67si/2 

67pi/2 

69si/2 

68p3/2 

2.84 

7.14 

-36000 

15.3 

3 

68si/2 

68p3/2 

70si/2 

69pi/2 

4.60 

7.34 

-23500 

13.1 

3 

69si/2 

69p3/2 

71si/2 

70pi/2 

-6.92 

7.80 

17600 

11.7 

1 

(M 

0 

00 

80pi/2 

CO 

00 

82pi/2 

5.20 

-14.9 

-76400 

15.7 

1 

81si/2 

81Pl/2 

84si/2 

83pi/2 

-2.92 

-15.7 

151000 

19.3 

1 

82si/2 

82pi/2 

85si/2 

84pi/2 

-10.3 

-16.5 

47400 

12.9 


TABLE IV. Dominant Forster resonances for Cs-Cs ns 1/2 states. The van der Waals interaction strength and Rc distance are 
given for the strongest eigenvector for the listed channel. 


channel 

Cs|7o) 

Gs 7a) 

Cs|76) 

Cs|7/3) 

5/2n (MHz) 

C 3 ,k (GHz /xm^) 

Uk, vdw (GHz /xm®) 

Rc (/im) 

4 

42si/2 

42p3/2 

43si/2 

42p3/2 

10.2 

-0.867 

-279 

5.49 

4 

43si/2 

43p3/2 

44si/2 

43p3/2 

-42.4 

-0.961 

82.2 

3.53 

2 

45si/2 

45pi/2 

47si/2 

46p3/2 

37.7 

1.28 

-86.3 

3.63 

3 

47si/2 

47p3/2 

49si/2 

48pi/2 

21.7 

1.41 

-184 

4.51 

1 

64si/2 

64pi/2 

68si/2 

67pi/2 

1.67 

-5.8 

-35700 

16.7 

2 

73si/2 

73pi/2 

76si/2 

75p3/2 

-1.88 

10.2 

110000 

19.7 

3 

76si/2 

76p3/2 

79si/2 

78pi/2 

7.59 

11.0 

-32100 

12.7 

3 

77si/2 

77p3/2 

0 

00 

79pi/2 

-3.91 

11.7 

69600 

16.2 

1 

81si/2 

81pi/2 

86si/2 

85pi/2 

3.86 

-15.7 

-114000 

17.6 


V. ANGULAR DEPENDENCE 


with the molecular axis R connecting atom a to atom b. 


The description so far has considered only the situation 
where the atoms are quantized along i which coincides 
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FIG. 3. (color online) Energy defects for interspecies cou¬ 
pling with nscs < nsRb, npnh = nsnh - 1, npcs = nscs- 
The different curves show the (jp,Rb,ip,Cs) channels in the 
sequence (3/2, 3/2), (1/2, 3/2), (3/2,1/2), (1/2,1/2) from top 
to bottom. 



FIG. 4. Energy defects for intraspecies coupling for the chan¬ 
nel nsns njp{n — l)jp for Rb-Rb (left) and Gs-Gs(right). 
The middle two series in each plot, which are almost identical, 
are for jp,jp = 1/2, 3/2 and 3/2,1/2. 



FIG. 5. (color online) Energy defects for Rb-Rb coupling. 
The different curves show the (jpa^jpb) channels. 
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The more general case of R at an angle 0 with respect to z 
is important for calculating interaction strengths in three 
dimensional ensembles. As we will see the near spherical 
symmetry of the interaction which is known for coupling 
of indistinguishable atomic nsi /2 states is substantially 
modified when we consider distinguishable atomic states. 

The angular dependence of the interaction is found by 
noting that when the molecular axis is rotated relative 
to a fixed laboratory frame the expansion coefficients Cki 
depend on the rotation angles, as described in [a , see 
also [ 2 ^. For the case of initial nsi /2 Rydberg states we 
have 


FIG. 6. (color online) Energy defects for Gs-Gs coupling. The 
different curves show the (jpa^jpb) channels. 


The angular behavior of the vanderWaals interaction 
in channel k is then found by generalizing Eq. (|3|) to 


^|V^),vdw(^) = ^ Uki,vdW 

i 

elk 1 


= -fk{0)- 


h5k 


(4) 


CkeiO) = 

1 /2 

with drn^,m^,{0) the reduced Wigner matrix for j = 1/2 
evaluated at angle 0. 


where fk{0) = \cke{0)\‘^Dke- 

The small R resonant interaction can be treated ana¬ 
lytically when a single channel is dominant. The interac¬ 
tion strength in channel k at angle 0 is given by Eq. (ED 
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with Dki replaced by fk{0) 


Uk{0) = ^ 


1-1 


h?5lR^ 


1/2 


(5) 


In the van der Waals limit we add the interaction from 
all channels. For s states there are two limiting cases 
of parallel and antiparallel spins. For the parallel initial 
state IV^) = 11 / 2 , 1 / 2 ) and 7 ^ 7 ^ 7 ^ we find 

Cki{0) = 0, 

Ck 2 { 0 ) = a/2cos( 6>/2) sin(6>/2), 

Ck3{S) = COS^(0/2), 

CkA{0) = s\v?{e/2). 


which results in 


/2 


/l = 


h = 


h = 


4 + 6sin^(0) 

9 ’ 

14 — 6 sin^(0) 
9 ’ 

22 + 6sin2(6») 


(6) 


For the antiparallel state = |l/2,—1/2) we find 


CfciW = l/\/ 2 , 
Ck 2 { 0 ) = cos( 6 »)/\/ 2 , 
Ck 3 { 0 ) = -sin( 6 »)/ 2 , 
Cki{9) = sin( 6 ')/ 2 . 


which results in 


/2 


/l = 


/s = 


/4 = 


8-6sin2(6») 

9 ’ 

10 + 6 sin^(0) 
9 ' 

26 — 6 sin^(0) 
9 ' 


(7) 


For other initial Zeeman states the angular functions fk 
will be different than those given here. Identical initial 
states with 7 ^ = 75 have fk twice as large as those in Eqs. 
(EEl). Although all channels have substantial anisotropy, 
the total interaction summed over channels, ^]^fk{^)i 
has no 0 dependence which implies that the interaction 
becomes isotropic in the limit of vanishing fine structure 
splitting between npi /2 — '^^ 3 / 2 - 

We proceed to illustrate these results with some 
examples. Consider two Cs atoms in the Rydberg 
state 81 si/ 2 |^a = 1 / 2 ), 81 <Si/ 2 |m 5 = 1 / 2 ) as a func¬ 
tion of 6. The interaction is dominated by the cou¬ 
pling |8l5i/2)|8l5i/2) ^ |81p)|80p) with all four fine 
structure channels contributing. Using C\^jbk = 
(-140.4,-237.7,-196.8,-463.0) GHz pwf for A: = 1,4 
we find 


ttCs-Cs 

^vdW 


3740. + 225. sin2(0) 

736 ((jHzj 


2.0 r 



°-° 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

0 (radians) 


FIG. 7. (color online) Angular variation of the van der Waals 
interaction for Rb-Rb, Cs-Cs, and Rb-Cs coupling at R = 
12.7 /xm. In each case the initial states have rria = rrib = 1/2. 


with R in /xm. For a pair of Rb atoms 
in the same 81 <Si /2 state we find C\^jbk — 
(-266.8,-362.1,-334.4,-499.1) GHz /xm^ and 




Rb-Rb 

vdW 


4840. + 92.5sin^((9) 
R6 


(GHz). 


The Gs-Gs and Rb-Rb interaction shifts are shown in Fig. 
[7|as a function of the molecular axis angle d. We see that 
both species have comparable interaction strengths that 
are weakly anisotropic with an angular variation of about 
2 % for Rb and 6 % for Gs. 

The situation can be markedly different for Rb-Gs. 
Taking ncs = 81 there is a Forster resonance with 
Rb at riRb = 77. The interaction is strongly domi¬ 
nated by a single channel k = 3. Using C\^jbk = 
(-729.0, 638.4,-69020., 345.2) GHz /xm^ for A: = 1,4 we 
find 


77Rb-Cs_ 105900.-45330. sin70) , 
l^vdW ~ ^6 (GHz). 

The interspecies interaction is stronger by about a fac¬ 
tor of 20, than for Rb-Rb or Gs-Gs, and is strongly 
anisotropic with a minimum at ^ = 7 r/ 2 . This is because 
of the dominance of the k = 3 channel. A different situa¬ 
tion arises for the Rb 84s 1/2 - Cs 89s 1/2 resonance. In this 
case Ci j^/Sk = (-136600., 1028., 2025., 683.9) GHz /xm^ 
for k = 1,4. The k = 1 channel is now dominant and 

7rRb-Cs _ 54300. + 92600. sin2(0) , 

l^vdW ~ (GHz). 

As seen in Fig. [7] the interaction now has a maximum 
at 0 = 7 r /2 and is strongly anisotropic. The Rb-Rb or 
Cs-Cs Forster resonances given in Tables innivi can also 
be anisotropic depending on which channels dominate. 
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VI. QUANTUM NONDEMOLITION STATE 

MEASUREMENTS WITH LOW CROSSTALK 
USING INTERSPECIES COUPLING 

One of the outstanding challenges of neutral atom ap¬ 
proaches to quantum computing is the requirement of 
qubit state measurements without loss, and with low 
crosstalk to proximal qubits. Such a capability is es¬ 
sential for implementation of quantum error correction. 
The most widely used approach to qubit measurements 
with neutral alkali atoms relies on imaging of fluorescence 
photons scattered from a cycling transition between one 
of the qubit states and the strong D2 resonance hne[2[. 
Due to a nonzero rate for spontaneous Raman transi¬ 
tions from the upper hyperfine manifold there is a limit 
to how many photons can be scattered, and imaged, with¬ 
out changing the quantum state. This problem is typi¬ 
cally solved by preceding a measurement with resonant 
“blow away” light that removes atoms in one of the hy¬ 
perfine states. The presence or absence of an atom is then 
measured with repumping light turned on, and a positive 
measurement result is used to infer that the atom was in 
the state that was not blown away. 

This method can indeed provide high fidelity state 
measurements but has several drawbacks. An atom is 
lost half the time on average, and must be reloaded and 
reinitialized for a computation to proceed. Atom reload¬ 
ing involves mechanical transport, and thus tends to be 
slow compared to gate and measurement operations. In 
addition, error correction would require that a single site 
in a qubit array can be reloaded, without disturbing prox¬ 
imal cm bits. While progress has been made towards this 
goal[2J-[^, much work remains to be done. 

Lossless quantum nondemolition (QND) measure¬ 
ments that leave the atom in one of the qubit states, 
or at least in a known Zeeman sublevel of the desired hy¬ 
perfine state, can be performed provided that the mea¬ 
surement is completed while scattering so few photons 
that the probability of a Raman transition is negligi¬ 
ble. This was first done for atoms strongly coupled to a 
cavity p7l - [^ . and was subsequently extended to atoms 
in free space [1, [ 30 .1^. 

Despite these advances, achieving useful QND state 
measurements in an array of neutral atom qubits re¬ 
mains an outstanding challenge due to the absorption 
of scattered photons by proximal atoms. Since the res¬ 
onant cross section for photon absorption is cr = 
and qubits in recent lattice experiments are spaced by 
d 5A|ii the probability of a scattered photon being 
absorbed is 77abs ^ cr/(47r(i^) ^ 0.0015. If the qubit mea¬ 
surement is performed with a moderately high numerical 
aperture collection lens of NA = 0.5 and the optical and 
detector efficiencies are 50% the probability of photon de¬ 
tection is T^det ^ 0.034 so that r^abs/^det ^0.04. This ra¬ 
tio implies that a state measurement based on detection 
of only a single photon would incur a ^ 4% probability 
of unwanted photon absorption at a neighboring qubit. 
This 4% error rate is too large to be efficiently handled 


by protocols for quantum error correction. 

One approach to solving the crosstalk problem is to 
protect nearby qubits in states that are dark to scattered 
resonant photons. This has been used effectively in ex¬ 
periments with trapped ions[32|. Such methods are in 
principle possible with neutral atoms, and an example 
using a single species is shown in Fig. [8] for ^^Rb. Sim¬ 
ilar ideas could also be implemented with other species. 
While the protection protocol can in principle solve the 
crosstalk problem it has the drawback of requiring both 
local and global operations, and is thus both complicated 
to implement and likely to be relatively slow. Neverthe¬ 
less this protocol points to an alternative approach using 
interspecies coupling. The Fig. [8] protocol suppresses 
crosstalk by placing all but the atom of interest in a dark 
state with respect to the probe light. Another way of sup¬ 
pressing crosstalk is to use one species for computational 
qubits and a second species for measurement qubits. Se¬ 
lective mapping of computational to measurement qubits 
allows us to probe the measurement qubits while keeping 
the computational qubits in a dark state with respect to 
the probe light, which is only resonant with the second 
species. 

This idea is made explicit using a two-species array 
as shown in Fig. [9l Our approach is analogous to 
the demonstrations of quantum logic spectroscopy [Tl| 
and entanglement [ 33 I with two ion species, and builds 
on earlier ideas of mapping single atoms to ensembles 
for fast readout [33 as well as the availability of asym¬ 
metric Rydberg interactions for creating multiparticle 
entanglement . The interspecies protocol requires 
fewer operations than in Fig. [51 and increases the useful 
photon rate per atom by a factor of four or more while 
eliminating crosstalk to other qubits. 

Consider the qubit array shown in Fig. |9l We as¬ 
sume this is a 2D array of 3D traps for Cs atoms as de¬ 
scribed in [il , and used for recent experiments with single 
qubit Q and two-qubit Q quantum gates. We will modify 
the array slightly by changing the wavelength of the trap 
light from 780 nm to 820 nm. This is still blue-detuned 
for Cs atoms which will be trapped at local minima of 
the optical intensity. The 820 nm light is red detuned 
for Rb atoms which will be trapped at local maxima of 
the intensity, forming a checkerboard pattern of alter¬ 
nating Cs and Rb atoms. The lattice period separating 
atoms of the same species will be d = 4 /im, and each 
Cs atom is surrounded by four Rb atoms at a distance 
of d = 4 /a/ 2 = 2.8 /im. The large wavelength separation 
between the Rb resonance lines at 780, 794 nm, the trap 
light at 820 nm, and the Cs resonance lines at 852, 894 
nm allows for independent loading, cooling, control, and 
measurement of the two species. 

Let us now choose an interspecies Forster resonance 
that gives strong Rb-Cs coupling and relatively weak Rb- 
Rb coupling. An example is shown in Fig. [10] for the 
Rb48s-Cs5l5 channel 3 resonance from Table [III Each Cs 
atom interacts with its nearest neighbor Rb atoms with 
a strength of t/RbCs/27r = 96.8 MHz at R = 2.8 /im and 
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global site specific global 

initial shelving measurement restore 



FIG. 8. (color online) Protection protocol for site specific measurements without crosstalk. Panel a) shows qubits in a 
superposition of hyperfine clock states of ®^Rb. In b) a global shelving operation is performed to map |a) onto |2,2) for all 
qubits. In c) a site specific mapping of \b) to |1,0) is performed. Resonance fluorescence is then generated using light of all 
polarizations coupling 5si/2|/ = 1) to 5p3/2|/^ = 0). The atom can only decay back to / = 1 and scattered photons are 
off-resonant with the shelved atoms in |2, 0). If scattered photons are detected the state is projected into j — \ and the 
measurement result was qubit state \b). The atom is then pumped into |1, 0) using cr+ and cr_ polarized light, after which it is 
rotated to |2,0) using a tt polarized microwave, or Raman light. If no scattered photons are detected the measurement result 
is qubit state \a) and the atom shelved in |2, 2) is rotated back to |2, 0) using a- — cr+ Raman light. Finally, in d) the globally 
protected states are restored back to \a) using Raman light. 


computational qubit measurement qubit 





o o o 





X=820 nmtrap array 


FIG. 9. (color online) Trap array with 820 nm light creates 
a checkerboard pattern of Gs computational qubits in blue 
detuned traps and Rb measurement qubits in red detuned 
traps. 


0 = 7r/2. We have assumed that the Cs and Rb atoms 



FIG. 10. (color online) Log-log plot of the interspecies cou¬ 
pling strength (upper curve) and Rb-Rb coupling strength 
(lower curve) for Rb48s and Gs51s states at ^ = 90 deg. The 
Rb-Gs curve is the full Forster interaction of Eq. ([2| for the 
single dominant channel 3 from Table HIl The Rb-Rb curve is 
the van der Waals interaction summed over all four channels 
for 48s48s o 47p48p. 


are excited to opposite m values so we use Eqs. (EEl) 
to calculate the interaction strength. In contrast two Rb 
atoms interact with a much smaller t/RbRb/27r = 2.3 MHz 
at = 4 jam. To measure the state of a Cs atom qubit 
we prepare all Rb atoms in the state |l)Rb = |2,2) by 
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optical pumping and then perform the sequence 

Cs : co|0) + ci|l) ^ co|0) + ici\bls) 

Rb : |I) ^ i|48s) —10) 

Cs : co|0) + ici\bls) -ir co|0) + ci|l). 

Provided the Rabi frequency of the Rb Rydberg exci¬ 
tation is small compared to t/RbCs we create the 
entangled state 

co|0)cs|0)Rb + ci|l) Cs |l)Rb- 

The overbar in the Rb kets denotes that this is a mul¬ 
tipart icle state of four Rb atoms. We then measure the 
hyperfine state of the Rb atoms. A detector click projects 
the Cs qubit into |1) and no click projects into |0). 

This approach has several advantageous features. 
Since each Cs atom is strongly coupled to four nearest 
neighbors the photon rate can be four times greater than 
for measurement of a single atom. This reduces the dif¬ 
ficulty of obtaining a hyperfine state measurement with¬ 
out suffering a Raman transition. Furthermore since the 
state of the Cs qubit is measured using fluorescence light 
at 780 nm which is far detuned from the Cs resonance 
lines, crosstalk to other Cs qubits will be negligible. After 
a measurement the Rb atoms can be rapidly repumped to 
the |l)Rb state in preparation for the next measurement. 
In addition to measurement of a single qubit the Rb 
atoms could also naturally be used as ancilla qubits for 
syndrome extraction in quantum error correcting codes. 

We proceed to estimate the measurement fidelity with 
realistic experimental parameters. When the Cs atom is 
in state |0) the transfer of the Rb atom between states 
will be affected by residual couplings t/RbRb to nearby 
atoms. This gives a transfer error for four atoms of 
E = 0.72I/^^j^^/|r2Rbp. The other dominant error is the 
imperfect blockade when the Cs qubit is in state |1). The 
first TT pulse on the Rb atom creates the state CRb|48<s) 
with 

If we set the Rb Rabi frequency such that t7RbCs/|^Rb| = 
^/3 then CRb = 0, and there is no state transfer, as de¬ 
sired. This result will be modified slightly by the pres¬ 
ence of more than one Rb atom, but an equivalent nulling 
condition will still exist. For the interaction strengths 
given above this condition is ORb/27r = 55.9 MHz and 
the transfer error is E = 0.72 x (2.3/55.9)^ = 0.0012. 
There is also a spontaneous emission error from the fi¬ 
nite lifetime Rydberg states. The Cs qubit is on average 
Rydberg excited for t = (l/2)(27r/f2Rb + 7r/f]cs)- The Rb 
atom is on average Rydberg excited for t = (l/2)7r/f2Rb. 
The room temperature lifetimes are[36| rRb48s = 58 /is 
and rcssis = 63 /rs. Taking ftcs = ^Rb we And Pse,Rb = 
7.7 X 10-^ and Pse,cs = 2.1 x 10“^. 

The largest error is the Rb state transfer at 0.0012. 
This small error occurs on average half the time when 


the Cs atom is in the |0) state and could be reduced even 
further by using a 25% larger lattice spacing which would 
increase the PubCs/PRbRb ratio by a factor of two. It is 
also likely that adiabatic or composite pulse sequences 
can be designed to minimize the sensitivity to small vari¬ 
ations in coupling strength}^. 

Finally we note that the use of two different species, 
combined with optical tweezers at a wavelength that only 
perturbs one species at a time, provides a means to move 
quantum information about in a larger array. This idea 
was developed for the case of Cs and Li atoms in Ref. 
[38| . In the cited work the entanglement of Cs and Li 
atoms was envisioned to occur via short range molec¬ 
ular interactions. The inter species Rydberg interaction 
described here can in principle be extended to Cs-Li, or 
other combinations, with the advantage that interactions 
can be performed at long range. 


VII. SUMMARY 

We have calculated the interspecies Forster interaction 
between Rb and Cs atoms, as well as Forster interac¬ 
tions for Rb-Rb and Cs-Cs where the participating atoms 
are excited to ns states with different principal quantum 
numbers. These interactions can be remarkably strong 
leading to van der Waals interaction strengths of several 
MHz at P = 20 jam for n < 90. The strong interactions 
are of interest for long range coupling between atoms of 
the same species which has already been demonstrated 
in Rb ensembles [^. 

We also propose to use the Rb-Cs interaction for loss¬ 
less and crosstalk free QND measurements. Needless to 
say the fidelity of this approach to measurements relies 
on having high fidelity Rydberg gates available. The cur¬ 
rent state of the art using the Rydberg blockade interac¬ 
tion, without post selection, uses a CNOT gate to create 
Bell states with a fidelity of 0. 730 . This is much lower 
than the intrinsic fidelity of the Rb-Cs mapping protocol 
which we estimate in Sec. l^to be ^ 0.001 with realistic 
experimental parameters. The two-qubit gate fidelity is 
therefore the largest roadblock for the protocol analyzed 
here. On the other hand, there is little interest in QND 
measurements of single atoms in a qubit array if high fi¬ 
delity gates are not also available. When a high fidelity 
Rydberg gate is demonstrated, the cross entanglement 
protocol described here may prove valuable for scaling 
up quantum information tasks with low cross talk. 
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Appendix A: Channel eigenvalues 

To find the angular factors il/c(ma,m^) for channel k 
and initial Zeeman states rria , we form the matrix of 
coefficients 
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/ 0 

0 



M = 

0 

0 

iy—3oc,-3p 

h5k 

lyjocjp 

0 


i /M-iad/3 

lyjoc ,i/3 

0 

Mk / 


(Al) 


The matrix has dimensions N x N with N = Nab + A^a/3 
and accounts for the coupling between states with the 
same value of m = rria + rub. The laser excited states 
are referred to as “initial” states and the dipole coupled 
Rydberg states as “target” states. The number of initial 
states is Nab = 1 + (ia Tib) — |^|- The number of target 
states Nai 3 is at most {2jo, + l)( 2 j /3 + 1 ), but may be 
less than that due to the requirement that rria mp = 
rria ^rrib. 


The nonzero off-diagonal entries are the dipole-dipole 
matrix elements 

^20 r^jo^mo, r^jprnf3 

q=-l 

(A 2 ) 

with Cs^k defined in Eq. O- The No,^ diagonals have 
value M/c- The eigenvalues and eigenvectors of M give 
the molecular energies of Rydberg excited atom pairs via 
a single interaction channel as a function of the atomic 
separation R with the quantization axis along R which 
points from atom a to atom b. 

When the j are half integers, which is the case for 
alkali atoms, and Nab = 1 the eigenvalues of M are of the 
following form. There are N — 2 degenerate eigenvalues 
U = hSk which have no R dependence and correspond 
to admixtures of a and /3 states. The remaining two 
eigenvalues are 


Uk± = 


hSk 


1± 1 + 4 


2^mc,=-jc, 2^mi3=-ji3 V^-^rna,mb ) 


1 / 2 - 


(A3) 


At large R the Uk- eigenvalue asymptotes to zero and therefore corresponds to [//c of Eq. ([ 2 ]) whereby we see that 

ia i/3 / \ 2 

A{m„,m,) = 6 £ ^ • (^4) 

ma=-jcy rni3=-ji3 \ q / 


When 7 a = 75 the eigenvalue is 2Dk. When Nab > 1 the is not possible to give compact expressions for [//., In 

eigenvectors are superpositions of \ma^mb) states and it these cases we extract the Dk from the calculated eigen¬ 

values by comparison with Eq. 0 . 















